Probabilities, Tensors and Qubits 



A. Yu. Vlasov 



26 April 2001 



o : 

\ Abstract 

^ ' In the paper is discussed complete probabilistic description of quan- 

turn systems with application to multiqubit quantum computations. In 

■ simplest case it is a set of probabilities of transitions to some fixed set of 
' states. The probabilities in the set may be represented linearly via coefB- 

\^ , cients of density matrix and it is very similar with description using mixed 

■ states, but also may give some alternative view on specific properties of 
quantum circuits due to possibility of direct comparison with classical 

, statistical paradigm. 

^ . 

^ ' 1 Introduction 

Let us consider n-dimensional Hilbert space Tin and set of N fixed unit vectors 
T— I ■ \va) G Tin denoted as S{va ■ a = 1,...,A^} or Sjv('Wn)- Then for quantum 
' system in arbitrary state \tp) G Hn there are defined N coefficients: 



Oh! 



The Pa is probability to find in state \va) due to measurement described by 
projector: 

^ ; n„ = n{va) = K)(va\. (2) 

O^' For mixed state with density matrix p it is possible to use instead of Eq. (IT]) : 



Pa = {Va\p\Va) =TT:{pna). (3) 



^ For pure state with p = the Eq. (g) coincides with Eq. (|l|). 

■ Due to Eq. probabifities of transitions Pa{p) are simply linear functions 
of density matrix and for properly chosen set of vectors Va using only the coef- 
ficients Pa may be equivalent to description by density matrices. For example 
in the paper is shown that for such a set any quantum channel can be described 
as linear transformation of iV-dimensional vector of probabilities, i.e. simply 
as NxN matrix. For n-qubit system it is possible to find a simple set with 
6" vectors, then the probabilities can be ordered in a tensor with n indexes, 
but only 4" components are linearly independent. It may be also described by 
4" special linearly independent parameters ordered in 4-dimensional real tensor 
with n indexes. A natural and symmetric description of any quantum circuits 
as linear transformations of these tensors is discussed below in this paper. 



2 General description 

Let us consider N probabilities Pa Eq. as formal vector p = {pi, . . . ,pn)- 
Let H{n) is space of Hermitian nxn matrices H G H{n), H^^ ~ H^^ and 
Hph{n) is "physical" subspace of nonnegative definite matrices with trace one, 

Hph{n) = {p : pG^f(n), Trp=l, {v\ p\v) > 0, V v}, (4) 
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then Sn via Eq. (||) produces Hnear maps Ls'. HpY^^n) K>q or Ls: H{n) 
. Here H{n) is hnear space described by real parameters and Hp\i{n) is 

subspace of physically possible density matrices p G i^ph(?^), Ls: p i— > P- 

Let us call the set Siv(H„) representative if exists right inverse of Lg, i.e., 

a linear operator LgiR^ — *■ H{n) with property: 

LsoLs=I, i.e., Ls(Lsff) = if, yHeH{n). (5) 

where symbol "o" is used for composition of operators (left one acts first). 

The space H{n) is described by real parameters and for minimal rep- 
resentative set with vectors S„2(7i„) the Ls is left and right inverse (or, 
simply, inverse: Ls = Lg ^) and together with Eq. (^) it is possible to write: 

LsoLs-I, i.e., Ls(Lsp)=p, VpeM"', (6) 

but if > A'^ the map Ls o Ls is not identity, it is projector on some n^- 
dimensional subspace V" C M^. The space V" is image of map Ls'. H{n) 

If representative set is not minimal, the Ls is not unique. There exists 
n'^{N — n^)-dimensional space of linear maps Lg- with property Ls"('^) ~ ^' 
Vv G V (here is zero Hermitian matrix) and so instead of Ls it is possible to 
use any other £.5 = Ls + Ls" . Due to such property results of alternative methods 
of construction of Ls may differ on if set S is not minimal. 

For physical space of density matrices Eq. (|^) due to condition Trp = 1, it 
is possible to express one of component using other with some affine relation. 
It is possible to introduce some affine minimal set with — 1 elements, but 
instead of linear equation Eq. (||) we will have affine equation: 

p = As(Lsp) + Ho, Vp e Hph{n). (7) 



Let us call a representative set Sat complete, if it can be expressed as union 
of few orthonormal bases of 7i„, i.e., for any Va G Ls it is possible to find other 
n — 1 vectors from Ls to form an orthogonal basis for 7i„. 

Let us call the complete set Sat almost perfect, if it can be expressed as 
disjoint union of few orthonormal bases of ?i„ and perfect, if for any Va £ Ls 
there is unique choice of n — 1 vectors from Ls to form the orthogonal basis for 
Ti-n- 

The definitions of complete, almost perfect and perfect sets are useful for 
consideration of elementary measurements with n possible outcomes. 

It is clear from definition that perfect set is also almost perfect. The example 
of almost perfect set that is not perfect, is one of main themes of this paper. It 
is complete set for m > 2 qubits introduced in Sec. ||. The set is almost perfect, 
but not perfect, due to counterexample in Sec. (page |ll|). It will be shown 
below (page ^) that almost perfect set may not be minimal. 

*** 

Let . . . , |n) is a basis of 7Y„ and pki are components of density matrix 
p = Pki\k){l\ and the map Ls: p 1— > p is defined as above by Eq. (H). Then 

the definition of Ls using Eq. corresponds to existence of n^N complex 
coefficients c^; (some of them may be zero and the choice is unique only for 
minimal set), where k,l — 1, . . . ,n and a ^ 1, . . . , N with property: 

N 

Pki — c'^iPa (there is no summation on k, I) . (8) 

Q = l 
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Let us introduce a basis on space C"^" of all complex matrices. Elements 
of the basis are Eki = |fc)(/|, i.e., matrices with only one unit in position (fc, I): 
{Eki)ij = SkiSij. Let us show, that a set is representative and Eq. is true 
for some coefficients c^; if and only if it is possible to express all r? elements of 
the basis as linear combinations of projectors 11^: 



N 

Elk = \l){k\ — ^ c^iHa (there is no summation on k,l) . (9) 

a=l 

First, Eq. (||) is direct sequence of Eq. because: 

N N N 

E ''kiPc = E Cfe; Tr(p n„) - Tr (p E ^ti^c.) - Tr(p \l){k\) = {k\ p \l) = pui 

a—l Q — 1 a— 1 

Now let us show that for representative set always exists some decomposi- 
tion Eq. (^). It maybe not obvious, the Eq. 1^ looks too general, because we 
considered basis of whole matrix algebra C"^" with complex (i.e. 2n^ real) 
parameters instead of space of Hermitian matrices H{n) described by only 
real parameters. 

The condition for representative set is existence of map Ls, such that com- 
position Ls o Ls is identity H{n) — > H(n) in diagram: 

H{n) ^MJ^ ^ H{n), Lg o Lj = I. (10) 

It is necessary and enough to satisfy the identity for whole basis of H{n) with 
elements. It is possible if > and in this case image of Ls must be some 

n^-dimensional subspace V" C M^. 

On the other hand, the Eq. (||) is equivalent with extension of the Ls, Ls and 

identity of Ls o Ls on whole space C"^": 

(^«xn ]l^^N }^ ^nxn^ L^ O L^ = L (11) 

But let us show that identity of Ls o Ls for basis elements H{n) described in 
Eq. ( p^ is enough for identity of the map for basis of C"^", i.e., the Eq. (11) 
is not more general. 

A basis of H{n) is n{n + l)/2 matrices 

H+ - [Eki + Eik)/2, k > I (12) 

(i.e. -tf^ = Ekk) and n{n — l)/2 matrices 

H- = i{Eki - Eik)/2, k>l. (13) 

But if some linear map C"^" (^nxn jg j^igjitity on these elements, it is also 
identity on linear combinations, but any matrix of basis E^i of C"^" may be 
expressed via and Hj^i : 

r H+~zH-, k>l 
Eki = < H+^, k = l (14) 

So Lg o Lg is identity and Eq. (||) is necessary and enough, it is simply 
expression for Lg written in basis Ei^- 

The following proposition summarizes the consideration above. 
Three properties are equivalent: 
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1. Set S(wct) is representative (Ls exists). 

2. Any complex matrix may be represented as linear combination of projec- 
tors r\{va) with complex coefficients — it is enough to show for any basis, 
for example Eki, see Eq. (^). 

3. Any Hermitian matrix may be represented as linear combination of pro- 
jectors r\{va) with real coefficients — it is enough to show for any basis, 
for example above. 

For minimal representative set the construction of Ls (item |^ in the propo- 
sition above) and both decompositions (items ^ and ^ are unique. 

Now it is possible to show that (almost) perfect set may not be minimal. 
For ?i-dimensional Hilbert space minimal set must have only elements. If it 
could be (almost) perfect set, it can be considered as disjoint union of n bases 
S{va(k,j)}i where k = 1, . . . , n is label of basis and j = 1, . . . , n is number of 
element in fc-th basis. But in this case we would have n decompositions of same 
matrix (unit) I = ^a(k,j)i ^^^d so it may be no more than — n + 1 

different linearly independent combinations of the projectors ^aik.j) instead of 
71^ and the set may not be representative. 

The complex basis Eki is "tr^s'-orthonormal" , i.e., orthonormal with re- 
spect to scalar product on space of matrices defined as 

{{A,B)),=TiiAB*). (15) 

Due to it there is simple relation between expression for the basis Eq. and an 



expansion Eq. (|^) of L for representative set. For Hermitian matrices Eq. (15) 
is close related with Eq. (||), because two different scalar products coincide: 

{{A,B))^TTiAB) = Tt{AB*) {{A,B)),. (16) 

The basis of Hermitian matrices defined by Eq. (|l2|) and Eq. (|l3|) is not tr^s- 
orthonormal, but it is possible to use standard Gram-Schmidt procedure of 
orthogonalization for given norm (cf. Ref. ||^). An example of tr^s-orthonormal 
Hermitian basis for qubits with dimension of Hilbert space is 2" will be discussed 
later. 



Other application of E^i is following theorem. 

Composition Theorem: Let 71^^ and Hm'^ are two Hilbert spaces with 
dimensions n and m and S]\[{va S Ti-n ^} and Sm{up € TiSi} are representative 
(complete, almost perfect) sets with N and M vectors. Then set of NM vectors 
'^NM{va®up G 'Hn^^'H.'^} is representative (complete, almost perfect) set for 
composite system. 

Let us prove the theorem for representative sets, because for complete sets 
it is trivial implication. For a proof it is enough to use condition Eq. (^ for 
n^n^ and n'^i^: 

N M 

\h){kA = Y.€Z^M. \h){k2\ = 5]cg,^n(K^), 

a=l 13=1 

(where |fci), \va) £ H^^'' and IZ2), |u/3) G H^^^) together with properties: 

n{v®u) = n(t;)®n(u), \iii2)(kik2\ = \ii)(ki\ ® |;2)(fc2|. 
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So for n^m^ elements of basis in space of all matrices C"™^""* = C"^"®C™^'" 
we have: 

N M 

\iii2){hk2\ = E E ^w^i.^i^c.mp). (17) 



The Eq. (17) corresponds to Eq. (|S|) for nm-dimensional Hilbert space 
Ti-nrn^ = Hn^ ^n!^ with representative set Snm{T^^^'^^) with NM vectors, if to 
use "compound indexes" hke \kik2\ ^ (fci — l)m + fc2 and \al3\ ^ (a — l)Af + /3: 



NM 



where cg^^ , , = ChC'i'J and = n(t;„®M^). 



A proof that tensor product of complete (almost perfect) sets is complete 
(almost perfect) is directly implied because tensor product of two bases is basis. 
Tensor product of perfect set is almost perfect, but not necessary perfect, it 



follows from counterexample in Sec. 5.2 



It is also possible to suggest proof without complexification, because tensor 
product of Hermitian matrices is also Hermitian, but it is not discussed here due 
to rather technical points. Such approach may be more clear from consideration 
with Pauli matrices in Sec. |^. 



3 Some applications of representative sets of 
projectors 

The Eq. (|^) shows that for representative set using of vector p is equivalent to 
description of quantum system by density matrix. 



Example 1: For any unitary operator U on 7i, |-!/;') — U\ 

p' = Up U* 



and 



(19) 



there exists NxN matrix Ajj, p' = Aup, there p' = p^/ (here Ajj does not 
depend on Really, the operator U induces a linear transformation U®U* 



Eq. ( |19[ ) on space of Hermitian matrices. Let us denote it as U 
Au = Ls o U*-'-"- o Ls as follows from diagram^]: 



Ls 



Ls 



Then 



(20) 



Let us consider construction of matrix Ajj more directly. Components of the 
vector p' are written as: 

p'„ =Tr(?7pC/*n„) -Tr(pt/*n„t/) (21) 

and because due to third item of proposition in Sec. ^ any unitary matrix can 
be represented as linear combination of V\a, we have 

N 

U*U^U = Y,AiUp, (22) 

n /3=1 

'^See also footnote (□) below on page M 
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where for given a, A\ , . . . , are real coefficients of decomposition of Hermitian 
matrix U*V\aU. The coefficients may be unique only \i N = and 11/3 are 
linearly independent, i.e., for minimal representative set. On the other hand, 
are components of desired matrix Au^ p' = Au'p, because 

N N N 

p'^ = Tr(p ^ Ain,) = 5] Tr(pn^) = ^ Aipp. (23) 

/3=1 /3=1 /3=1 

If set is not minimal, the decomposition is not unique, for example it is 
possible to use only n? linearly independent projectors between N > n? 
and already this choice is not unique. There is a problem with calculation 
of coefficients A^ even for minimal set, but it is convenient for set of matrix 
orthogonal in respect of some norm. The complex extension with matrices Eki 
orthogonal in norm Eq. (15) was discussed already and for system (n = 2"*) 
of 771-qubits basis of Hermitian matrix obtained from 4'" different products of 
Pauli matrices orthogonal in both norms Eq. (16) and Eq. ( |l5|) is used below 
Eq. (H). 

Let us suggest, that we choose some basis of Hermitian matrices and 
express all matrices 

Ha = 51 h^HK. (24) 

K=l 

The coefficients are unique and if basis Hk is tr^iB-orthonormal 

T,{HjHk) = 5jk (25) 
then the coefficients may be expressed as 

h^K = Tr(n„ffj^) (26) 
Let us introduce new parameters: 

PK = THpHk). (27) 
Then it is possible to express all N > probabilities pa using the parameters: 

2 2 2 

n n n 

p^ = Tr(p n„) = Tr(p ^ H^Hk) = ^ Tr(pi?K) = ^^P^' (^8) 

K=l K=l K=X 

If set is not minimal, coefficients /i^ of inverse transformation: 

N 
a = l 

and 

AT 

PK = J2 ~^KPa (30) 

a=l 

are not unique. 

The parameters pK can be considered as a n^-dimensional vector p, but now 
transformation p' — Aup always unique and can be calculated directly, let 

p'K^TriUpU*HK) = TripU*HKU), (31) 
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it is possible to use decomposition 



U*HkU^Y.^kHj, (32) 
j=i 

and so 

p'k^T.^kPJ, (33) 
J=i 

where matrix Au is expressed as: 

{Au)i = Ty{HjU*HkU). (34) 

Initial matrix Ajj can be expressed from Au using linear maps and h'^ 
between pa and px- If set is not minimal, Ajj depends on h'^ and is not unique. 
It is unique only restriction of An on ri^-dimensional linear subspace V discussed 
above. It is possible also directly work with pK, because all probabilities can be 
expressed using the parameters and expressions Eq. (|2^) where are always 
unique. The px can be considered as "coordinates" on V. 

The idea to write transformations of state of quantum system in terms of 
probabilities (weights) was also suggested in relation with some other problems 
in ||, ^ , but authors found linear transformation only for n = 2 and work 
only with pure states.^ 



Example 2: The representation may be even more convenient for general 
quantum channel [0 : 

p'^Y^^kpV^, Y,VkV^=l. (35) 

fc k 

In this case linear transformation^ V: p' p is defined by some set {Vk} and 
the set is even not unique, but because the diagram Eq. ( po[) is valid for any 
linear map H{n) H(n), any quantum channel can be expressed using only 
one matrix = X^fe^Vt, p' = ^vP- For minimal set matrix Ay is unique, 
overwise unique only restriction of A\[ on n^-dimensional linear subspace V. 

Applications of parameters p here is also justified. It is directly followed 
from consideration of previous example. A matrix of map A^: p i— > p' is also 
always unique. 

^Initial version of present paper (May 2000) appears as positive answer on question if the 
maps like An may be linear for n > 2. It should be mentioned also, that for pure state due 
to polynomial relations between coefficients of density matrix like pkiPmn = Pknpml it may 
be possible to express some probabilities in minimal set from the others by some nonlinear 
functions and so use nonlinear transformations with lesser amount of vectors Va instead of L 
and Ajj, but it is not discussed here. On the other hand, even for linear Ajj the existence 
of inverse map Lj for Lj is too strong requirement. If to look only for the linear map Ajj, 
it is enough to consider condition that diagram Eq. ( po| ) is commutative. In this case it is 
reasonable to consider irreducible representations of tensor product U{n) X U*{n). It can be 
shown, that such product is sum of two irreps with dimensions are — 1 and 1. The second 
one is simply scalar representation corresponding to constancy of density matrix trace and 
so splitting off the one-dimensional irrep may be explained without applications of theory_i3f 
group representations. Due to such situation, instead of commutativity of diagram Eq. ( poj ) 
here was used straightforward, but less general suggestion about existence of Ls, especially 
because the map is anyway necessary for other applications discussed in this paper. 

^Sometimes it is called "superoperator." 
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Example 3: Any expressions for probabilities can be found using only ele- 
ments of vector p (cf. Ref. 1^). For any vector v ^ Ti there is linear map 
defined as 

{v\p\v) ^ By^ = {v\Ls^\v) (36) 

{By depends only on f, not Let p^ corresponds to some pure state t/i, 

p^ = Ls(|V')(V'l) by definition of hsv via Eq. (0), Eq. (|), then Eq. ^ can 
be rewritten as |(f | V')!^ = B^Pff,. It is probability of transition. Let p^ = 
Ls(|w)(w|), then because \{v \ = Tr{\v){v\\ij){i;\) = Tr (Ls(p.)Ls(p^)) is 
some bilinear map X ^ R, due to general theorem of linear algebra, it 
is scalar product (Pu,p^)s — \ {v \ ip)\^ defined by some matrix Gs'- 

N 

(P,q)s= E Gfpaq0. (37) 

a, (3=1 

So it is possible to define some metric Gs on with property that probability 
of transition for pure states is scalar product in this metric. 

Example 4: ft is possible to calculate average value of any Hermitian operator 
X for some observable: X''' = {^\X\^P) = Tr(Ls(pv.)^) (cf. Refs. §, § 
§8.1.1]) for pure state. For mixed state it is possible to use formula X'' = 
Ti{pX) — Tr(Ls(pp)X) and if to represent a Hermitian operator X as some 
vector px = Ls(X) g M.^ , then average value of operator X for some mixed 
state p is expressed via the same metric Gs Eq. (^7|) introduced in example 
above 

(px,Pp)s. (38) 



4 Example of representative and complete sets 

Let us consider example with Hilbert space Tin and minimal representative set 
with N ^ n? elements. The -n? vectors are expressed using n basis vectors \a) 
as following three families of vectors: 
n basis vectors itself: = |a), 

vectors: = + |/3)), a < /3, 

vectors: K^) = ^(|a) + a < p. 

The set is known also due to application for construction of POVM pl| . 
The weight vector p with N — n? real components is composed as set of 
numbers from three different families: 

Pa = (<^l P l"^) = Paa 

Pip = \{{a\ + {P\)p{\a) + \f3)) = \{paa + Ppp + {pa[3 + Ppod) 

PaP = U{a\~-i{l3\)p{\a)+i\l3)) = y^Paa+ Pl3p+i{Pal3- Ppa)) 
So paa = Pa and {pap + P/3q)/2 = pj^^ - 5 (p^ + P^p) , ^{PaP + P/3q)/2 = 

Pap ~ ^(Pa ~^ Pp) for a</3 and it is possible to write for arbitrary matrix: 

{p^ ; (y — j3 

plp-^plp-^^ip^+pp); a<P (39) 
Pla+^Pla-'-^iPl+P})-^ 
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or for Hermitian matrix it is convenient also to use expressions with real 
parameters for real and imaginary components of elements pap (here a < [3): 



Pact — Vai 

The Eq. (^9|) is direct construction of linear map L discussed above and so 
r? vectors are representative set. 

But the set is not complete. It is possible to add two other families of vectors 
for completion: 

vectors: = - a < /3, 

"^T^^ vectors: |u2'«) = 45(|a) - «|/3)), a<(3. 



5 Complete sets for quantum circuits 
5.1 One qubit 

The complete set described in the previous example has 2n^ — n elements and for 
case of qubit n — 2 with basis |0), |1) G H2 it is especially simple and symmetric 
perfect set with six elements: 

|0^) ^ \vl) = |0), \V) ^ \vl) = 

|0-) ^ \vf,) = ^(|0) + ^ = ;^(|0) - |1)), 

|0^) ^ K2) = ;^(|0) |F) ^ K%) = ^(|0) ~ z|l)). 

It is simple to check that they are eigenvectors of Pauli matrices with eigen- 
values ±1: 

fj^'lB'') = AbIB'") (no summation), Ag = (-l)^ (40) 

where 8 = 0,1 and p G {x,y,z}. The Eq. (^ ) can be simply checked using 
expression for projectors: 

= n(0'^) = |0^)(0^| = i(I + As a"), (41) 

= AB(2n^ - 1) = 2n(^ - 1 = I - 2n5' = - n^", n^ + n^ = i, (42) 

then 

a^|B^) = AB(2n^ - I)|0^) = AB(2n^|0^) - |0^)) = Ab(2|0^) - |0^)) = Ab|0^). 

Because operator of spin can be expressed^ as ^cr^ the Eq. ( ^ ) shows that 
\0^) and |1'') correspond to spin +i or — ^ respectively for measurements in 
three orthogonal directions p G {x,y,z}. 

Six equations for probabilities p§, . . . ,pf are: 

= (BHpIB'^) =Tr(pn^) =Tr(pi(I + ABa^)) = i(l + Ab Tr(pa'')). (43) 

Let us now use together with p £ {x, y, z} numerical indexes p — 1,2,3 and 
also p,iy = 0, . . . ,3: 

(T — IL, (7 — a , a = T" , a = CF . 

''in Planck's units 
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It is a basis of space of Hermitian matrices: 

3 

jy = ^a^a'', HeH{2), a^eW. (44) 

Because of a property of the Pauli matrices: 

Tvia^a'') = 26^, (45) 

they form the tr^^-orthogonal basis and it is simple to find the coefficients 
in Eq. (H): 

= iTr(i?(T''). (46) 

For density matrix it is always qq = ^, but here is used all four indexes 
due to application to composition below. Let us denote — ~ p'^ and 
po = = 1, \JfjL = 1, 2, 3. Then from Eq. (||) follows: 

p,^Ty[pa-), z/ = 0,...,3. (47) 

Up to insignificant multiplier i, pi, correspond to parameters an just defined 
in Eq. ( ^ ) or pK introduced earlier in Eq. (27). The density matrix can be 
expressed as 

3 

p=lY,p^a\ (48) 

due to Eq. (^ij) and probabilities can be expressed as 

= 5(Po + Aap^), Ai= 1,2,3, = 0,1, po = I (49) 

due to Eq. (H). 

A simple transformation property of parameters (here p, — 1,2,3) for 
unitary 1-qubit gate (cf. example 1 in Sec. ||and Q) is related with 2-1 homo- 
morphism S'0(3) and SU{2) groups — for any matrix U € SU{2) there is 3x3 
matrix O G 5*0(3) with property: 

3 

= ^O^^d''. (50) 
1^=1 

For a map Ajj -.p^ ^ P'^^ we have: 



p'^ = Tr{UpU*an = Tr(pC/V^C/) = Tr(p^O;>"^) = ^0;Jp.. 

So for = (pi,p2,P3) it is possible simply to write: 

p' = O^ip, Ou e 50(3). (51) 

The Eq. ( [5l|) shows that one-qubit gates corresponds to 3D rotations of vector 
(pi,P2,P3) M and parameter po is not changing. 

Here is used mainly a model with spin-half systems, but an analogy of the 
parameters Eq. (^), Stokes parameters was also considered in optical models 
of qubits pi. 
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5.2 Two qubits 

Let us apply the composition theorem for system of two qubits. Here we have 
complete set of 36 = 6^ vectors 185*^92^) with ^i,fi2 — 1,2,3. The set is almost 
perfect^ it can be considered as disjoint union of nine bases with four elements: 
{|0''iO'^^), \0^'n^'^), ll'^iO^^), for nine pairs In example with 

spin systems, it corresponds to nine possible combinations of measurements of 
two spins {S , S fj^rj) with four possible combinations of results, (±i,±i), i.e. 

The set is not perfect, for example for vector |0^0^) the choice of basis 
{|0^0^), |0^1^), |1^0^), |1^1^)} is not unique, it is simple to check that vectors 
{|0^0^), |0^1^), |1^0^), |1^1^)} are also orthogonal. The example shows that ten- 
sor product of two perfect sets may be almost perfect, but not perfect. 

It is possible to consider 36 (linearly dependent) probabilities: 

Ki^^ =Tr(pn^^^®n^^^), /ii,M2 = 1,2,3, Qi,B2 = 0,l (52) 
described by 16 parameters: 

p,,,,^Tr{pa'''®a'''), i/i, = 0, 1, 2, 3. (53) 
The Hermitian basis with 16 elements ^a'^^^a'^^ is trAs-orthonormal (see 



Sec. 5.3 below) and it is possible to write 

3 

1/1, 1/2=0 

The expressions of p via p can be simply calculated 

<^^2 = Tr(pn^^^®n^,^)-iTr(p(a° + Aoia''0®(^° + AB2'T^^) 

iPfJ-lO + AB2_P0/i2 + Agj Ab2P^i/J2)- (^5) 

where poo = Tr p = 1 . 

On the other hand, to express p via p it is necessary to use few kind of 
formulas (here /ii,/J.2 ^ 0): 

p^,^, = Tr(p(n(;^-nf)®(n(;^-nf)) = K^(;^-p(;^f-pf[;^+pff, 
P0M2 - Tr(p (n^^ + nf )®(n^- - nf )) = p^^f;^ - p^o'^ + p^i'^o' - f> 
p,,o = Tr(p (nj;^ - nf )®(n^^ + nf )) = p^^r +p^^f - - pf f, 
1 - POO - Tr (p (n^^ + nf )® (nf + nf ) ) - pf f + pf f + pf f + pf f, (56) 

there only for nine parameters p^i/^a with /ii, /X2 ^ it is unique decomposition, 
for six parameters with one zero index p^^^o or Poai2 there are three different 
expressions and the last equation shows nine different ways to decompose unit. 
It is possible to introduce six probabilities for description of first subsystem, 

p^'^' EETr(pn^(g)I), /i = 1,2,3, B = 0,1 (57) 

and six for second one, 

p^'^^ = Tr(pl®n^), M=l,2,3, B = 0,1. (58) 

The equations correspond to measurement of only first or second spin re- 
spectively for models with two spin-half systems. 
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Each such probabihty can be expressed as sum of two probabihties described 
by Eq. (52) in three different ways (corresponding to free parameter /zi or ^,2)'- 

P^'^'^=K^^+P^f, Va*i,M2 = 1,2,3. (59) 

It is also possible to introduce similar description with parameters p, but in 
this case it is only change of notation: 

=Tr(pa''®I) =ii,o, = (0, ) 1, 2, 3, (60) 

=Tr(pI®a'0 =Po., (0,)1,2,3. (61) 

(where p"^^ = poo = Tr p = 1, Vj). It is also possible to use the parameters 



Eq. ( pOD and Eq. (|61|) to express probabilities Eq. ( pTf ) and Eq. (|58| ) respectively: 

= i(l + (-l)^pW), J = 1,2, M = 1,2,3, 8 = 0,1. (62) 

Let us consider situation without entanglement, when: 

P = Pi®P2- (63) 

(cf. \^) = \'4'i)®\'4'2) for pure states) 
Using identities: 

{A®B) {C®D) = {AC)®{BD) (64) 

together with 

Tr(A(8B) = Tr(A) Tr(S) (65) 
we have for any matrices Ai, A2'. 

Tr((pi®p2) {Ai®A2)) = Tr((pi Ai)®{p2 A2)) = Tr(pi Ai) Tr(p2 ^2). (66) 

So if systems are not entangled, then due to Eq. ( p3| ) and Eq. (^) it is 
simple to prove expression for "independent probabilities" . Due to definitions 
Eq. (ph, Eq. (|7|), and Eq. (pi), it is possible to write: 



P'ilZ = Tr(pi®p2 n^^^^n^^^) = Tr(pi H^^^ Tr(p2 ^Z)^ 
Psl = Tr(pi®p2 n^^^^I) = Tr(pi H^^^) Tr(p2) = Tr(pi H^^), 
pZ = Tr(pi®p2 I^n^) = Tr(pi) Tr(p2 H^) = Tr(p2 U^), 



and, finally, 

p','X^p's'/'^p',:^'^ (67) 

Similar expression 

P.:.. =P«P('^ (68) 

can be proved for definitions Eq. (|5^), Eq. (|60|), and Eq. (^. 
5.3 Many qubits 

For composite system with m qubits density matrix is p G iJph(2™) and it is 
possible to use 4™ different tensor products of Pauli matrices as a basis: 

3 

J2 a^,...,.r^<j'''(E}---(g)cr'''", H e H{2"'), a^,...^,„ £ M. (69) 



Let us denote cr'^i - '^'" = a"^® ■ ■ ■ ^a""^ . Due to Eq. (|65|) it is possible to 
write: 

Tr(^.i....™^m...p™) = 2™<5,,^, •••<5,„^„, (70) 
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i.e., the Hermitian basis 2 '"cr'^i - '^™ ig tr^s-orthonormal. An expression for 
coefficients a,j-^...u„^ directly follows from the Eq. (|71 



= 2-™Tr(iJCT'^i-''™). (71) 



So for density matrix aoo...o — 2^™ Tip = 2^™. 

Using composition of m perfect qubit bases with six vectors discussed above 
we can produce basis with 6™ components. It is almost perfect because can be 
considered as disjoint union of 3™ bases with 2™ elements: |B^^ • • ■ 0^") marked 
by 3™ different sets (/ii, . . . ,^m)- It is similar with description of case m = 2 
above and the counterexample provided there is enough to show that the set is 
not perfect for any m > 1. 

Let us denote flg^^ = Ha^^ig) • • • ^Hg". When 6™ probabihties introduced 
by complete (and almost perfect) set S6™(Hf"") can be represented as: 

<:::b,: = Tr(pn^^^:::^;j, Aife = 1,2,3, 9,. = 0,1, k^i,...,m, (72) 

described by 4™ parameters: 

=Tr(pa'^i-'^'"), i.fc = 0,l,2,3, fc = l,...,m, (73) 

where po...o = Tr p = 1. 

In example with spin systems, Pb^"'b^ correspond to 3™ possible combina- 
tions of measurements of m spins (5^^ , . . . , >S'p„ ) with 2™ possible combinations 
of results {^Xbi, ■ ■ ■ , ^As^). 

It is possible to introduce six probabilities for each subsystem, similarly with 
Eq. (It]) and Eq. (||), 

pM(fe) = rj,^(^j®fe-i^(-|^^jj»m-fc>)^ ^==1,2,3, = 0,1, /c = l,...,m (74) 
described by three parameters similarly with Eq. (|60| ) and Eq. (|6l|), 

pW=Tr(pI«'=-i®a''®I«'"-'=)=Po...o.o...o, k^l,...,m, (75) 

fc-l m-k 

(where = po...o = Trp = 1, Vj). The Eq. ( |6^ ) again can be used for 



calculation of probabilities Eq. (74), 

i(l + (-l)apW), k = l,...,m p= 1,2,3, = 0,1, (76) 

but analogue of Eq. ( |59| ) is too difficult, because it would contain 2™^-'^ terms 
and may be expressed in 3"'""'^ different ways. It is yet another example of 
usefulness of parameters p. 

Here is also can be considered systems without entanglement, 

m 

p = Pl® ■ • ■ = (^Pfc. (77) 

fc=l 

For such states is also can be used analogue of Eq. ( |67[ ) 

m 

P^l::Z ^l[Ps:^'^ Mfc = 1,2,3, 0fc = O,l, k=l,...,m, (78) 

k=l 

and Eq. ( |68| ) 

m 

= n^S^ ;^fc = 0,l,2,3, fc = l,...,m. (79) 

fe=i 
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So, systems without entanglement can be described by 6m probabilities 
Eq. d?!) instead of 6™ Eq. (||) or 3m parameters Eq. (|7|) instead of 4" - 1 
parameters Eq. For systems with entanglement such simple idea does 



not work and 6"* probabilities of Eq. ( 72 ) and 4™ parameters Eq. (|73|) can be 



considered as some tensors, Pg = Pg^\\\gl^ and Pi, = p^^...u^. 

Transformations of the tensors due to action of local quantum gates on initial 
state is also local. If some quantum gate acts only on a few (Z = 1, 2, . . .; / ^ m) 
qubits and can be described as 2' x 2' complex matrix, then transformation of 
tensor can be described by 4' x 4' real matrix (and 6' x 6' matrix for pg ). 

It is clear after rewriting of expressions with trace like Eq. (|3^) already used 
above, i.e, 

pIi....,„ = Tr({/pC/V"i'-'^'") = TripWa-^'-'-^U) (80) 

and so if U acts only on few indexes between z^i, . . . , v^, then only the indexes 
present in matrix Ajj of the transformation p' — Aup or, formally. 



E 



Aii::Zlp,....,^ (81) 



Really it follows directly form expressions for coefficients Eq. ( p^ ) that can be 
rewritten as 

^tl'"!^!,") = Tr(cr'^i-'^"[/V^i-^'"[/). (82) 
For example for two-qubit gate acting on some indexes j and k instead of 



Eq. (81) for Ajj with 4™ components it is enough to use only 4 — 16 com- 



ponents and write: 

3 

fvi...Vj...Vk---Vm / , (UjV,,)t^Vi...Hj...^lk..Mm- \"^J 

The similar approach is true for Ajj and pg . 

Such probabilistic description of quantum circuits is not more complex, than 
usual description with pure states, and for description of quantum channels with 
mixed states it is even more simple, because here is not necessary to use some 
special decomposition for "superoperators" , they are also may be expressed by 
one matrix with same size as for "usual" operators. 

Here again was rather used model with m spin-half systems, but parameters 
similar with was also used as "multiple beam Stokes parameters" in rather 
optical framework jl^ . 



6 Conclusion and bibliographical notes 

The set of problems discussed in present paper is part of more general task of 
complete description of quantum systems using some set of observables. The 
task had a long history and very wide range of different approaches and it is not 
possible here to review that, but together with references used in main body of 
the paper it should be mentioned some alternative approaches |@, |[ || together 
with references therein. Present paper itself appears initially in rather specific 
circumstances of discussions related with macroscopic realization of quantum 
logic |2| 1^, ^, ^ with later development in more general framework related with 
quantum computations. 
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